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Decay channels with attractive branching ratios, or interesting physics, are recovered by substi- 
tuting "missing" particles (7, f, vr", etc) with combined geometric and kinematic constraints. The 
"Sliding Vertex" method is shown in this part-I, for reconstructing strongly boosted decays - at 
the LHC. 

PACS numbers: 13.20.He,14.40.Nd 
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Precision B-Physics in today's experiments [H relies 
heavily on exclusively reconstructed decay chains. There 
are however decay modes that are not readily exclusive- 
reconstructable, but very attractive from the view point 
of the physics, or that of a large branching ratio - im- 
paired in reconstruction usually by missing particles. 
Methods for constraining and recovering information in 
missing particle events have been explored before with 
kinematic fitting [3, Q and in a variety of other contexts, 
a good collection thereof being [3, Q . 

Consider for instance the B° — > A'+ decay with 
the subsequent K+ tt^ (tt^) decay. This mode 

is simple and attractive for measuring 7ckm ; however its 
branching ratio is rather fair. If we consider the related 
mode, with 7r° in the final state, the branching ratio is 
ca. 3.5 times larger, a real feat. Evidently, the 7r° can be 
reconstructed in the E-calorimeter, albeit with less reso- 
lution than for tracking. In principle the reconstruction 
with the tt'^ is very attractive, but seemingly somewhat 
impractical. 

There is however enough information in the detector to 
reconstruct the decay without the 7r°: first there are the 
4- momentum conservation laws in the Bvtx and Dyxx, 
and second, the Byxx must lie on the 7r+ track. Figure[T] 




FIG. 1: Topology of ^ D^ K'^ decay with the subsequent 
DJ K+ K" TV- {tt°). The tracks are shown in the absence 
of magnetic field - such as for LHCb, although, given the small 
track bending over the vertexing region, a solution in the form 
of an Oiii helical correction exists also for ATLAS. 



shows the topology of the event (in the absence of mag- 
netic field - such as for LHCb, although, given the small 
track bending over the vertexing region, a solution in 
the form of an Om helical correction exists also for AT- 
LAS). The evident question arises, whether such meth- 
ods reject not only background, but also "sister physics" 
(in the case of the above mode the significantly more 
abundant B^ —^ T)J tt"*" decay, with the tt"*" mistaken by 
the PID system as a K"*"), which play a more dominant 
role. Such is the case of decays coming from B" which to 
a large degree shadow (ca. 3:1) the B° decays through 
similar mass ratios between mother and daughters. From 
the start it should be stated that this is the most diffi- 
cult part and that possible solutions lie in the avenues 
of tracking/vertexing resolutions, of (semi)-leading (bi)- 
particle effects in the aforementioned decays and of dif- 
ferent branching ratios (sometimes smaller) for B^ for 
the exact event topology as B^. In fact B° contaminates 
B° events more heavily with different topologies (where 
1 track is lost, or which have 1-2 extra tt^'s). 

Kinematic Constraints - are the 4-momentum con- 
servation laws in the Bvtx and Dvtx • For the Byxx the 
useful part is: 
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and for the Dyxx^ 
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where "0" is the missing neutral, "v" the sum of visible 
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particles in the Dyxx and the rest of notations are self- 
evident. 

Geometric Constraints - are supplying the missing, 
3'''^ equation to the above set. Due to momentum conser- 
vation in the Byxx "D", "tt" and "q" (see figure 1) lie in 
the same plane: nu ^ An^r + M^q, where n are unit vec- 
tors and A, fj, constants. The cosines from the kinematic 
relations are: cosOott = ^Z5 • nn and cosOuv = nu ■ nv, 
respectively: 

A = [{vq)cos9D7r - {qn)cos9Dv] 

A* = [ - {vtt)cos9dtt + cosOdv] (3) 



where (ab) 



fi}, and A = [vq) — (w7r)(g7r). The 



"closure" equation is \nD\ = 1 = A^ + //^ + 2{q7r)\fi. 

In terms of the geometric constraints the kinematic 
section condenses to E"^ — = (Xpo)^ + (fJ-Po)^ + 
2{qT:)(XpD){p-PD) which can be solved in favor of Ed as 
a second order equation. The two fold ambiguity result- 
ing therefrom is lifted through a (2D) pointback to IP 
criterion for the (after vertex determination). 

At this point all kinematic quantities are knowufs']. 

Vertex Determination - once E]j is known the 3D 
vectors of all particles are known. The ByTX b is the 
locus that takes the best-shot at the: 



We coded our own MXV4 [y] namespace that addressed 
these issues and performed in speed. 

One special mention is with respect to the tt track 
matrix. This is not immediate two-fold: writing the a,^ 
and secondly, the tube is not cylindrical, rather elliptical 
in cross-section. This is due to to the fact that the track 
is determined by (approximately) circular errors in the 
VELO planes of LHCb (prototype addresse) which stand 
vertically. (This of course remains to be solved exactly 
for each detector where it is applied.) The tt track tube is 
thus an infinitely long ellipsoid with unequal cross-section 
major axes: 
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where a.. 



and 



With this, the tt track contribution in the Bvtx calcu- 
lation reduces significantly (and this needs to be so im- 
plemented in the code in order to avoid singularities). 



1. IP in the ng direction: b — Xb^b = IP i c/p 

2. DvTX in the no direction: b — Xonn = -DvTxicrri 

3. 7r+ track Icttt tube: b — XT^fiT^ = tott ± ctd 
Mathematically this means - using the Ui error matrices: 

( b-VB - Xbtib I crjp \b-VB - Xbtib ) + 
{ b-VD- Xdud I ctd^ I b-VD~ XofiD ) + 

{ b -foTT ~ X^fijr I Cr^^ I & - rOTT - XjrnTr ) = TOm(4) 

By differentiating to find the minimum, the "sliding" 
along each direction is: 



Xi ^ { Hi \ a^^ \ b - Vi ) / { fii \ \ni ) 
The vertex solution is then: 
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The math in itself is straightforward, however com- 
puter implementation proved to be somewhat of a hassle. 
Code that implements vectors and matrices is rather slow 
due to multiple inheritances, unnecessary functions, etc. 
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FIG. 2: The xy- (top) and z-residuals (bottom) for the recon- 
structed BvTX with respect to the Monte Carlo position. All 
dimensions are in /im. 

We tested the code on PYTHIA simulated data 
(smeared to give invariant mass resolutions as those ex- 
pected in the LHCb detector). 

The xy- and z-residuals for the reconstructed Byxx 
with respect to the Monte Carlo position are shown in 
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FIG. 3: Momentum of tt" coming from (blue) vs. that 
from D~. The spectrum is harder for the latter, as the other 
two components do not necessarily each contain a valence 
quark of the decaying particle, leading to a more equitable 
energy distribution among daughter particles. This is not 
the case for the two kaons comming from a T)~ decay, which 
together take a greater share of the decay energy (leading- 
particle effect 01 • 

figure m With an rms of 130 /im and a central part 
of 30 /im, the z-rcsolution is very good compared to 
typical expectations on the order of 200-300 /im. The 
xyrms is 11 /im with a central core of 4 /im, an excellent 
value. As mentioned, it may be possible to reconstruct 
the mode, but then there are further problems associ- 
ated with background suppression, in particular that of 
similar-physics. In principle soft/QCD backgrounds are 
kinematically very far off and eliminable. The events that 
are a concern are those from Bj], insomuch not those of 
the same topology, rather those of different topology and 
1-2 lost tracks, and/or mis-ID tracks. 

Leading particle effects - figure [3] shows the tt^ mo- 
mentum (yellow) from — + tt^ tt^ K+ of the decay B° 
— > D~ TT"*", where tt^ is mis-ID'ed as a kaon. Charge 
does not help in distinguishing, as both B and B can be 
present, mass differences are small, hence the only factor 
that can still play some role would be one due to dy- 
namics. Since the CM decay energy is 0.7-0.8 GeV some 
leading-particle 0] effect is visible for the kaon-pair in 



a Dj decay (each holding a valence quark of the Dj). 
In blue is the tt~ momentum coming from Dj, evidently 
softer (less available energy, as most is concentrated by 
the kaon pair). 

Applications - the method aims evidently at physics 
analyses, however other useful events can also be recon- 
structed: ones that have neutrinos (used in tagging), 
or ones with many pions in one vertex (for instance B[J 
D*^ 7r+ TT^ 7r+ where the Bvtx has 4 pions). The 
purity of such pion samples make them ideal candidates 
for studies of tt — > K mis-ID rates. 
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